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ABSTRACT 

Context. The corona is highly dynamic and shows transient events on various scales in space and time. Most of these 
features are related to changes in the magnetic field structure or impulsive heating caused by the conversion of magnetic 
' to thermal energy. 

, Aims. We investigate the processes that lead to the formation, ejection and fall of a confined plasma ejection that was 

^ i observed in a numerical experiment of the solar corona. By quantifying physical parameters such as mass, velocity, and 



:3 
<N 



orientation of the plasma ejection relative to the magnetic field, we provide a description of the nature of this particular 
plasma ejection. 

Methods. The time-dependent three-dimensional magnetohydrodynamic (3D MHD) equations are solved in a box ex- 
tending from the chromosphere, which serves as a reservoir for mass and energy, to the lower corona. The plasma is 
heated by currents that are induced through field line braiding as a consequence of photospheric motions included in 
the model. Spectra of optically thin emission lines in the extreme ultraviolet range are synthesized, and magnetic field 
lines are traced over time. We determine the trajectory of the plasma ejection and identify anomalies in the profiles of 
the plasma parameters. 

Results. Following strong heating just above the chromosphere, the pressure rapidly increases, leading to a hydrody- 
^ ' namic explosion above the upper chromosphere in the low transition region. The explosion drives the plasma, which 

^. needs to follow the magnetic field lines. The ejection is then moving more or less ballistically along the loop-like field 

I . lines and eventually drops down onto the surface of the Sun. The speed of the ejection is in the range of the sound 

O ' speed, well below the Alfven velocity. 

^ ' Conclusions. The plasma ejection observed in a numerical experiment of the solar corona is basically a hydrody namic 

, phenomenon, whereas the rise of the heating rate is of magnetic nature. The granular motions in the photosphere lead 

. (by chance) to a strong braiding of the magnetic field lines at the location of the explosion that in turn is causing strong 

currents which are dissipated. Future studies need to determine if this process is a ubiquitous phenomenon on the Sun 
on small scales. Data from the Atmospheric Imaging Assembly on the Solar Dynamics Observatory (AIA/SDO) might 
, provide the relevant information. 

^ ! 
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OS 

• 1. Introduction often non-trivial and might be attributed to either waves 

or mass flows. For example, early on, intensity varia- 

O The solar corona is highly structured and very dynamic. ^-^^^ ^^^^ associated with the evaporation and conden- 

^ This ranges from large disruptive events such as flares and ^^^-^^ of plasma and the m otion of the condensation re- 

T-H coronal mass ejections (CMEs) down to small-scale explo- g^^^ (lAntiochos fc Sturrockl . flQTSi) . Outward-propagating 

^ sive events m the transition region between the chromo- features of fan-like structures have ofte n been interpreted 

. ^ sphere and the corona. Because the coronal dynamics are a ^s (slow) magneto-acoustic waves, (e.g., iDe Moortel et ID . 

direct signature of the processes governing the energy bal- [2002I) 
ance in the upper atmosphere, i.e., coronal heating, radia- 
tive losses and heat conduction, the investigation of coronal 

transients on all temporal and spatial scales is a key for also The investigation of cool plasma of several 10.000 K as 

understanding the nature of the corona of the Sun and other seen in the Ly-a line above the limb indicates that small 

stars. We present results from a three-dimensional magne- brightenings are moving along the mag netic field lines, 

tohydrodynamic (3D MHD) numerical experiment of the These have been analyzed in depth by ISchriiverl (|200l[ ) 

solar corona in which an explosion-like event drives a con- using data from TRACE (T ransition Region and Coronal 

fined plasma ejection ballistically through the corona along Explorer: Handy et al.Lll999[). Th e movie provided with the 

the magnetic field lines. online version of ISchriiverl ( 200 if ) clearly shows three types 



Imaging instruments provide information on the tem- of bright cool structures moving in the corona: (1) features 

poral and spatial variation of the coronal emission, appearing up in the corona that are falling down, (2) bright- 

Especially during transient events, their interpretation is enings moving upward that subsequently disappear, and (3) 

small brightenings initially moving upward before they fall 

Present address: International Space Science Institute, back again. These types of brightenings have been inter- 

Hallerstrasse 6, CH-3012 Bern, Switzerland preted in different ways. 
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Type (1) features can be interpreted as condensations in 
the corona caused by a loss of thermal equili brium. These 
features hav e been investigated in depth by iMiiller et alJ 
(|2QQ3l |2QQ4) through ID loop models and compared to ob- 
servations (|Muller et al.Ll2QQ5f ). These models show that if 
the heating is concentrated in the lower atmosphere and 
is not sufficient to balance the radiative cooling at the 
loop apex, a runaway process will be initiated, because the 
cooler plasma can cool more efficiently and through this 
rapidly forms a condensation. Eventually, this condensation 
slides down the field lines with speeds of up to lOOkm/s. 
iDe Groof et al.l (|2QQ4l ) observed that features moving down- 
ward are indeed plasma flows rather than magneto-acoustic 
waves. 

The second type of cool brightenings that move up- 
ward and subsequently disappear might be associated 
with the hi gh-speed transition region line upflows as dis- 
cussed by iMcIntosh fc De Pontieul (|2QQ9[ ). The authors 
conclude that these upfiows might be visible as prop- 
agating disturbances in the extreme ultraviolet. They 
may or may not be the small-scale relatives o f exp losive 
events and spicu l es. [Mcintosh fc De Pontieul (|2009) and 
iDe Pontieu et all (|2QQ9[ ) found these upfiows through a 
blueward asymmetry and line broadening of transition re- 
gion and coronal lines in active regions that have been 
report ed be f ore for tran sition lines in the quiet Sun net- 
work (|Peterl . |2QQQ[ l2QQlf ). The cool brightenings eventually 
heat up and disappear from the observations of the cool 
plasma (as in Ly-a), leaving a footprint in the asymmetry 
of lines form ed at higher coronal temperatures above 10^ K 
(jHara et al.L [2008; Peter[ j2010[ ). 

The type (3) cool brightenings that move up and fall 
back might be interpreted as plasma that is propelled up- 
ward, follows the magnetic field lines and then falls back 
onto the Sun without being heated up to coronal tempera- 
tures as the type (2) brightenings are. We found evidence 
for this type (3) process in our numerical experiments of 
the solar corona. The driving of the magnetic field in the 
photosphere through horizontal convective (granular) mo- 
tions results in a high energy input over a short time caused 
by Ohmic dissipation in the transition region. An explosion 
event of this kind would be pushing cool material upward, 
which then follows the magnetic field lines up to the apex 
and slides down on the other side of the magnetic loop. 

The high energy deposition that eventually leads to 
the plasma ejection comes about self-consistently when (by 
chance) the granular motions in the photosphere lead to a 
magnetic configuration with a high rate of Ohmic dissipa- 
tion in a small region in the simulation box. 

By analyzing the data of the numerical experiment, we 
will investigate the formation, rise and fall of the plasma 
ejection in the 3D MHD model and quantify its physical 
parameters, such as mass, velocity and orientation relative 
to the magnetic field. In addition, we synthesize emission 
lines that are observable with current instrumentation. 

We emphasize that this study is not investigat- 
ing spicule-like phenomena. Related high-speed upfiows 
that have recently been observed from the chromo- 
sphere through t he t ransition region to the corona 
(|De Pontieu et all . [201l[ ) most likely play a significant role 
for the supply of hot plasma into the corona. Likewise, ob- 
servations have s hown strong upfiows at the edges of ac- 
tive regions (e.g.. [Sakao et al.' . . [2007[ : iDoschek et al.l [2QQ1 . 
Observational evidence exists that these upfiows appear 



as asymmetries in the line pro f iles (e . g . , iDe Pontieu et al.[ . 
I2QQ9I : IMcIntosh fc De Pontieul . l2QQ9D_ or as second compo- 
nents in the spectral profiles (e.g.. iHara et al.L 2008; Peterl 
|2010[ ). However, these high-speed upfiow structures seem to 
be relatively stable. Probably, the transient plasma ejection 
we discuss here is closer related to miniature c oronal mass 
ejections (mini-CMEs) as described by, e.g., Ilnnes et aD 
(2009, 2010.). 

We briefly explain the setup of the numerical ex- 
priment in Sect.[2l A qualitative description of the ob- 
served phenomenon is provided in Sect.[3l and the on- 
set of the eruption is investigated in detail in Sect.lH 
The dynamic evolution of the ejection along its tra- 
jectory is discussed in Sect.O and in Sect. [6] we com- 
ment on the visibility of this phenomenon with respect 
to the Atmospheric Imaging Assembly onboard the Solar 
Dynamics Observatory (AIA/SDO). Finally, the results are 
summarized and discussed in Sect. [71 

2. Setup of the numerical experiment 

We analyze an ejection of cool material into the corona in 
a 3D MHD simulation of the solar corona above a small 
active region. The s etup of the 3D model is identical to 
the one describe d inlBingert fc Peterl (2011). We used the 
Pencil code (Brandenburg fc Doblerl . [2002) to solve the 3D 
MHD equations. The heat conduction term parallel to the 
magnetic field and the radiative losses were included in the 
energy equation to obtain a realistic description of the coro- 
nal pressure and thus of the emerging synthesized emis- 
sion from the transition region and corona. The corona 
is heated through Ohmic dissipation of currents that are 
induced by the photospheric footpoint motions shuffling 
around the magnetic field. The computational domain en- 
compasses 50 Mm X 50 Mm horizontally and 30 Mm verti- 
cally. The initial magnetic field at the lower boundary is 
taken from an observed active region magnetogram. After 
approximately 30 min of simulation time, the solution is in- 
dependent of the initial condition, and the magnetic field 
in the corona is computed self-consistently based on the 
MHD equations. After this time period, the corona has 
reached a dynamic equilibrium, where the overall appear- 
ance and spatially averaged quantities do not change sig- 
nificantly, while the situation is still dynamic with strong 
flows and changes in the plasma properties (cf. movie at- 
tached to Fig.lA.l]). For more details, we refer the reader to 
iBingert fc Pet 61^20 111 ). 

The model does not contain the proper physics of the 
chromosphere, especially (non-LTE) radiative transfer is 
not included in the calculations. However, because the 
phenomenon we describe emerges from above the chro- 
mosphere, this is not expected to be a severe issue. At 
the top of the chromosphere, the energy input is typically 
10~^ Wm~^, while for the plasma ejection described here, 
the heating rate just above the chromosphere rises to val- 
ues a factor of 100 higher than that just before the ejection. 
Therefore, we do not expect the inclusion of non-LTE radia- 
tive transfer to have a major effect here. Additional energy 
transfer, e.g., through radiation, would certainly slightly 
mitigate the effect, but would not fully compensate the 
strong increase in heating rate (see also Sect. 14. 3p . 

From the results of the 3D MHD model, we calcu- 
lated the emission in the transition region and corona. 
Emission line spectra of various transition region and coro- 
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Fig. 1. Evolution of the plasma ejection over 30 minutes (top to bottom, times noted in right panels). The left and middle 
panels show the (normalized) synthesized emission in O VI (1032 A) and Mgx (625 A) integrated horizontally through the 
computational box. This corresponds to what would be observed at the solar limb. The right panels show the currents, 
j^, i.e., the Ohmic heating rate, also integrated horizontally through the box. is divid ed by the average trend with 
height to account for the roughly exponential decrease of the heating rate with height (cf. iBingert fc Peterl . l201l[ ). The 
arrows point to the same location in the maps in each row. The red line in the bottom row indicates the trajectory of 
the plasma eruption that was deduced from the position of its center-of-gravity (see section |42]) . The red box in the 
j ^-images shows the approximate formation region of the plasma eruption; the dotted lines indicate its extension in the 
x-direction. See also movie in Appendix O Fig. lA.li 



nal li nes are synt hesized us ing the CHIA NTI atomic data 
base (|Dere et al.l, |l99 7; Lan di et al.l , I2006D . This procedure 
follows IZacharias et al.l (|200ll201lD . The O VI (1032 A) line 
was chosen to represent cool transition region material, and 
the Mg X (625 A) line to represent hot coronal plasma (see 



Tabled]). The emission of these lines is roughly compara- 
ble to the 304 A and 171 A bands observed with AIA/SDO 
(cf. Sect. [6]), whic h were derived using the temperature re- 
sponse functions (|Boerner et al.L l201l[ ) as provided by the 
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Table 1. Emission lines and AIA/SDO bands synthesized 
from 3D MHD model. 



line 


wavelength 


formation temperature 




[A] 


log(T/[K]) 


Ovi 


1032 


5.5 


Mgx 


625 


6.0 


AIA / He I 


304 


4.9 


AIA / Fe ix/x 


171 


5.9 



Note 1: For O VI and Mgx, the line formation tem- 
perature is given, i.e., the temperature where G{T) oc 
(nion/riei) T~"^/^ exp [— /izy/(/csT)] reaches its maximum value 
(assuming ionization equilibrium). 

Note 2: For the AIA bands, the main contributing lines are given 
along with the center wavelength and the formation tempera- 
ture of the plasma dominating the emission in the respective 
band. 

instrument team in Solar Soft Q The plasma is assumed to 
be optically t hin an d in a state of ionization equilibrium 
(Peter et al., '2006^). The constant elemental abundances 
are taken from Mazotta et al. (1998). 

The specific choice of the O VI and Mg x lines is not 
crucial. The temporal evolution of the ejection looks quite 
similar in Civ and Nevill, respectively. The synthesized 
Mgx images are similar to those synthesized for the 195 A 
channel of AIA, the Nevill images are similar to those of 
the 171 A AIA band. In this study, the difference between 
these four coronal images is not significant, and therefore 
the validity of the images displayed in Figs.[TJ [9] and lA.ll 
does not depend on the choice of the respective line. Also, 
the small differences in appearance between images of C IV, 
Ovi and AIA 304 A are not significant for the illustrative 
purpose of these figures. Actually, O Vl is more similar to 
the AIA 304 A band than Civ. The AIA 304 A band has 
some contribution from high temperatures, which origi- 
nates from hot emission lines in the bandpass. Likewise, 
O VI shows some contribution from high temperatures that 
is caused by di-electronic recombination. C IV, being a Li- 
like ion, also shows this effect, but only to some extent. The 
quantitative analysis in Sects.|4]and[5]is not influenced by 
the choice of the lines or bands, of course. 

3. Phenomenology 

We are focusing on a time sequence of about 70 minutes 
duration in the numerical experiment. All points of time 
discussed here refer to time t=0 at the beginning of this 
sequence. The actual simulation was started long before 
t=0, so that any impact of the initial condition is excluded 
during the investigated time sequence. 

The ejection lifts off at about t=40min and lasts for 
about llmin. It is observed in the form of a "bubble" of 
plasma that rises and falls back down again along a curved 
trajectory. Snapshots of the emissivities integrated along 
the horizontal ^/-direction of the simulation box are pre- 
sented in Fig.[T] for different time steps of the simulation 
run. In addition, the corresponding Ohmic heating rates 
(oc j^, the current density squared) are shown. For details 
on the displayed images and a movie of the temporal evo- 
lution, see Appendix [Al and Eqs. lA.ll and IA.2p . 

^ http:/ /www. lmsal.com/solarsoft/ 




Fig. 2. Coronal box as seen in Ovi (1032 A) synthe- 
sized from 3D MHD model (emission increases from blue 
over violet to red). Overplotted are the magnetic field 
lines that cover the trajectory of the ejection (snapshot at 
45min). Visualiza tion using VAPOR (iClvne fc Rastl . [2005[ 
IClvne et all . 120071) . 



Before the plasma eruption, an increased heating rate 
is observed just above the chromosphere in the region out- 
lined by the red box and arrow a in Fig.[I] at t=38 min 
(best seen in the movie with Fig. lA.l) indicated by an ar- 
row). At t=41.5min, the plasma ejection lifts off (Fig. (TJ 
second row), and the cool ejection is clearly visible in the 
O VI line formed at transition region temperatures of about 
3-10^ K (arrow h). The cool ejection, however, is hardly 
visible in the Mgx coronal line formed at approximately 
10^ K (arrow c). In this line, the plasma ejection merely 
appears as a darkish shade because of its low temperature 
compared to the surroundings. At t=47min, the ejection 
reaches its maximum height and starts descending. When 
falling down, the ejection compresses the magnetic field at 
its front. As a consequence, currents increase in front of 
the ejection (arrow /), plasma is heated, and the front is 
brightening (arrows e; again, best visible in the movie). 
The magnetic field compression caused by the ejection pro- 
duces a heating front that is traveling ahead of the actual 
plasma ejection. Thus, the plasma is heated at low heights 
just above the chromosphere and consequently shines bright 
in the coronal line (arrow h). At this location and time 
(t=52min), the emission in the transition region line is not 
yet brightening (arrow g)^ but it does when the ejection ac- 
tually hits the chromosphere (as can be seen in the movie). 
There is a time delay between the brightening of the coro- 
nal line and the transition region line (at the location of 
arrows g and h) of about 1 min. In the movie, some of the 
plasma can also be seen to bounce up again in the rear part 
of the ejection because of the increased pressure when the 
ejection hits the chromosphere. This observation is simi- 
lar to the bouncing of condens ations falling down on the 
solar surface, as was found by iMiiller et all (|2003L I2004D . 
A few minutes later, the plasma ejection has completely 
disappeared (Fig.dJ bottom row). 

The correlation of the plasma ejection and the magnetic 
field is depicted in Fig.[2l which provides a 3D view of the 
computational domain as it would be seen in the emission 
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of the O VI line at time t=45 min, i.e., close (in time) to the 
middle row of Fig.[Tl 

The ejection, shown together with magnetic field lines 
in the vicinity of its trajectory in Fig.[2l has almost reached 
its maximum altitude and is moving nearly parallel to the 
magnetic field structure (cf. Sect. 15. 3p . As can also be seen 
in Fig. [21 the magnetic field is expanding with height, which 
causes the varying diameter of the ejection (cf. Fig.[T] and 
movie provided in AppendixlX]) . Fig. [2] also shows that the 
magnetic field is not helical. 

4. Onset of the plasma ejection 

In this section, the processes that trigger and drive the 
plasma ejection will be investigated. To study the onset of 
the plasma eruption in terms of time and space, peculiar- 
ities in the plasma properties must be looked for, mainly 
with regard to the current density, which is an indicator for 
plasma heating that can cause an eruption (Sect. l4?T]) . 

When looking at the computational domain from the 
side, the ejection becomes clearly visible after it has 
emerged above the ambient transition region. This is 
the case at a height of approximately 10.5 Mm, between 
t=39min and 40 min. The shape of the ejection appears 
tailed when ascending or descending and ellipsoidal when 
it is close to its maximum height. The ejection returns to 
the initial height at t=51.5min on the other side of the 
loop-shaped magnetic field lines. Subsequently, it slowly 
dissolves into the bright plasma and has completely dis- 
appeared after t=56 min. 

To determine the trajectory of the ejection, the plasma 
that is part of the ejection must be defined. It proved to 
be useful to consider ejected plasma in the volume where 
the emissivity of Ovi is observed to exceed 10~^ W m~^. 
Thus, the surface of the ejection can be defined and its 
center of gravity localized. Tracing the center of gravity as 
a function of time allows us to determine the trajectory of 
the ejection. This procedure works best for the time periods 
after the ejection has detached from the ambient transition 
region, during its rising phase and before it touches the 
ambient transition region again when it falls back onto the 
solar surface. The analysis of the plasma parameters along 
this trajectory is presented in Sect. 14. 21 

4.1. Heating rate and pressure in the formation region 

Fig. [3] visualizes the evolution of both the heating rate and 
the pressure in a small volume indicated by the red box 
in Fig.[T] (the volume extension is identical for both the x 
and the y direction) chosen to capture the formation site of 
the ejection. Before the eruption, the Ohminc heating rate, 
^Ohm=Mo^j^, continuously increases by approximately two 
orders of magnitude within 20 minutes (Fig.[3l left) causing 
a temperature increase in the formation region. At the same 
time, the density increases through chromospheric evapo- 
ration. These occurances lead to a considerable rise in pres- 
sure (Fig.S]) that eventually results in the observed explo- 
sion at t^39 min. 

The rapid rise of the Ohmic dissipation, and thus the 
pressure, is limited to the small volume indicated by the 
red box in Fig. [T] extending 2 Mm x 2 Mm horizontally and 
1 Mm vertically). In general, the heating of the corona in the 
applied model is caused by the braiding of magnetic field 
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Fig. 4. Propagation of the pressure disturbance in the early 
phase of the ejection. During this interval, the ejection still 
rises almost vertically. The different curves show the (peak) 
pressure as a function of height in a vertical column indi- 
cated by the dotted hues in Fig. [1] (bottom right panel). As 
time progresses, the location of the peak pressure moves 
upward. The vertical extension of the red boxes in Fig. [T] is 
indicated by the vertical dashed-dotted lines. 



lines (|Parkerl[T972lll988[ ). Numeric al 3D MHD models for 
this process have been developed bv lGudiksen fc Nordlundl 
(2002, 2005a.bi) for the first time. In these models, horizon- 
tal motions in the photosphere shuffle around the magnetic 
field which results in the braiding of magnetic field lines. 
Thus, currents are induced that are subsequently dissipated 
in the upper atmosphere. 

The photospheric driv ing by the horizontal m otions in 
our 3D MHD model (see iBingert fc Peterl . l201l[ ) leads by 
chance, not by prescription, to a strong increase of the heat- 
ing rate at the spot that will later become the origin of the 
plasma ejection. This rapid increase of the pressure is com- 
parable to an explosion that subsequently drives plasma 
away from the explosion site. 

According to the realistic driver used to describe the 
photospheric motions, it appears plausible that a strong 
increase of the current density just above the chromosphere 
can also happen on the Sun, e.g., in the periphery of an 
active region or pore. 

4.2. Pressure along the trajectory of the ejection 

The trajectory of the ejection must be defined before we 
investigate the ejection properties along the trajectory. We 
determined the volume of the ejected plasma as described 
above (i.e., where the emissivity of Ovi exceeds 10~^ W 
m~^) and calculated its center of mass for each snapshot 
of the numerical model every 30 s. Thus, the path of the 
ejection can be tracked from t=40 min to 51.5 min, at which 
time the entire plasma ejection is visible and distinguishable 
from the background. For time steps prior to t=40 min, the 
trajectory is obtained by a linear downward extrapolation 
(basically vertical). 

For the early onset of the ejection before t=39 min, the 
vertical column containing the formation site of the ejec- 
tion was examined in detail just below the spot where the 
ejection becomes visible. This column has a quadratic cross 
section (2 Mm x 2 Mm) and is indicated in Fig. [T] (bottom 
right) by the vertical dotted lines. The (peak) pressure in 
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time [min] 

Fig. 3. Temporal evolution of plasma properties in the source region of the plasma ejection. The panels show (a) the 
Ohmic heating rate i^ohm=Mo^j^7 (b) the work done by the magnetic field and (c) the work done by the pressure gradient. 
Positive values in panels (b) and (c) imply that work is performed on the plasma, i.e., the kinetic energy of the plasma 
is enhanced. All values represent averages in a box spanning 2 Mm x 2 Mm horizontally and 1 Mm in height (cf. red box 
in right panels of Fig.[T]). 



the vertical column is shown as a function of height for 
different time steps in Fig.lH The dashed-dotted vertical 
lines in Fig. |4] indicate the vertical extension of the red box 
shown in Fig[Tl 

At t=32.5 min, a small kink appears in the pressure pro- 
file at a height of approximately 5 Mm, indicating a pressure 
anomaly. Remaining more or less at the same location, it 
develops into a local pressure peak within a few minutes. 
This marks the phase of temperature increase and evapora- 
tion of the plasma in response to the increased heating rate. 
It is followed by the explosion, and, as a result, the pressure 
peak starts moving upward, as can be seen on snapshots at 
times t=37min to 39 min. The part of the pressure dis- 
turbance moving downward is effectively dissipated by the 
dense plasma of the low atmospheric layers and is therefore 
not visible. 

Typically, a disturbance that travels upward into re- 
gions of lower density steepens and eventually forms a 
shock. However, dissipation can counteract this tendency 
to form a shock, and thus, the shape of the perturbation 
will be invariant on its way upward in the simulation box. 
This phenomenon of a soliton-like structure is indeed what 
is observed here. The pressure disturbance appears in the 
form of a soliton-like structure that keeps its amplitude and 
half width, leaves the top of the red box at t=38.5 min and 
continues to travel upward. 

The subsequent evolution of the pressure of the ejec- 
tion on its way upward is shown in Fig.[5l As the maximum 
pressure of the ejection is found roughly at its center of 
gravity, these values match the ones shown in Fig.|4] (for 
better comparability, the actual height in the simulation 
box is shown rather than the arc length along the trajec- 
tory). Thus, the lowermost location of the ejection as shown 
in Fig. [5] at t=38.5min exactly corresponds to the upper- 
most location of the ejection shown in Fig. SI which also 
occurs at t=38.5min. This shows that the plasma ejection 
continues to travel upward, still in a soliton-like fashion, 
but now the dispersion slowly becomes more important: 
The shape flattens as the ejection expands in the course of 
the flight. 
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height [Mm] 

Fig. 5. Similar to Fig. SI but for the peak pressure of the 
ejection along its trajectory (shown as a dashed line in the 
lower panels of Fig.[TJ see also red line in Fig. [8] for a 3D im- 
pression). For better comparison with Fig. SI the pressure is 
plotted as a function of height, viz. altitude, along the tra- 
jectory (and not as a function of arc length along its path). 
The first time step (39 min, leftmost peak) corresponds to 
the last time step shown in Fig. SI 



4.3. Energy considerations for the onset of the ejection 

In the previous section (|4.2p we outlined a scenario in which 
the increasing heating rate leads to an increase in pressure. 
The resulting pressure gradient drives the ejection — sim- 
ilar to the situation of an explosion. We will confirm this 
scenario by investigating the energy budget. 

Fig. SI shows that the pressure p increases by 
0.3 Nm~^ within r=5min (the pressure starts increasing 
at t=32.5min and reaches its peak at about 37.5 min). For 
an ideal monoatomic gas, the internal energy density e is 
given through e=3/2p. Thus, the rate of change of the in- 
ternal energy density is about e/r = 1.5mWm~^. This 
nicely corresponds to the Ohmic heating rate i^ohm of 1.5 
- 2mWm~^ during that time (see Fig.[3l left panel). Thus, 
the pressure increase is indeed caused by the increased 
Ohmic heating rate of the plasma. The heating rate of some 
2 mW m~^ is about a factor of 100 higher compared to nor- 
mal regions at or just above the top of the chromosphere. 
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The acceleration during the early phase of the ejection 
can be derived from Fig. HI At times t = (38, 38.5, 39) min 
the ejection reaches heights of about z = (5.4, 5.7, 6.7) Mm 
when traced by the peak pressure. 

From this, one can calculate the velocity to be about 
V ~ 25kms~^ and the acceleration to be a ~ 800kms~^. 
The density in this region is p ^ 5 - 10~^^ kgm~^. Thus, the 
rate of change of the kinetic energy density of the ejection 
is about dt{\pv'^) = pva ^ O.lmWm"^. This number cor- 
responds well to the amount of work done by the pressure 
gradient as shown in the right panel of Fig. [3] and gives a 
clear indication that the work done by the pressure gradient 
is pushing the plasma upward. 

We can therefore conclude that the ejection is driven 
by the pressure gradient. The middle panel of Fig. [3] shows 
the work performed on the plasma by the magnetic field, 
or more precisely by the Lorentz force. This is about a fac- 
tor of 10 weaker than the work performed by the pressure 
gradient, and it would not be strong enough to drive the 
ejection. Therefore, the Lorentz force does not play a sig- 
nificant role for the onset of the plasma ejection. When the 
ejection leaves the source region just before time t = 40 min, 
V'{jxB) becomes negative, i.e., the plasma performs work 
on the magnetic field by pushing aside the magnetic field 
lines. 

Finally, we note that the kinetic energy of the plasma 
ejection is about 10^^ J or 10^^ erg (assuming 10^ kg for the 
average mass and 50kms~^ for the average velocity of the 
ejection). This is about a factor of ten higher than th e typ- 
ical energy content of a nanofiare (e.g., iParkeri 1 19881 ). 

5. Dynamics and evolution of the plasma ejection 

5.1. Mass balance and chromospheric evaporation 

After defining the volume of the plasma ejection, the total 
mass of the plasma ejection can be tracked. After a total 
of about 0.5-10^ kg has been ejected following the initial 
explosion, the gas bubble acquires additional mass of about 
10^ kg on its way along the magnetic field (see Fig. [6]). The 
total mass of the coronal part of the computational box, 
i.e., the mass of the plasma at temperatures above 10^ K, 
is found to be about 5-10^ kg. Therefore, the cool plasma 
ejected up into the corona is a significant fraction of the 
total coronal mass. Still, it is only a tiny fraction of the 
entire mass in the computational domain owing to the high 
density of the photospheric and chromospheric plasma. 

The (linear) increase of the mass of the ejection that 
flies through the corona as depicted in Fig. [6] cannot sim- 
ply be accounted for by an accumulation of plasma in the 
corona before the lift-off of the ejection. Before and dur- 
ing the ejection, the (curved) tube along the trajectory of 
the ejection is filled by evaporated chromospheric material. 
The evaporation is caused by the increased heating rate 
before and during the plasma ejection. Thus, without the 
occurence of the eruptive ejection it is likely that a coronal 
loop (with increased density) would have formed through 
the long-known process of chromospheric evaporation fol- 
lowing increased heating. After the explosion- type event, 
the ejection collects all material, which prevents the coro- 
nal loop from forming. The reason for this is that it sweeps 
away the evaporated material and transports it back to the 
chromosphere on the other side, where it is thermalized 
once the ejection hits the surface again. 




^ 0.5 



C/2 

I o.o F - 

40 42 44 46 48 50 52 
time [min] 

Fig. 6. Temporal evolution of the mass of the plasma ejec- 
tion. 




40 42 44 46 48 50 52 40 42 44 46 48 50 52 

time [min] time [min] 



Fig. 7. Left: Temporal evolution of the vertical velocity; 
overplotted in gray are linear fits in the time intervals 
t=40 — 45 min and t=46.5 — 50.5 min. Right: Temporal evo- 
lution of kinetic plasma- /3 at the center of gravity of the 
plasma ejection (see Eq.[T]in Sect. 15. 3p . 



5.2. Dynamics 

The vertical velocity Vz of the plasma ejection is shown 
as a function of time in the left panel of Fig. [71 The de- 
pendency is characterized by the vertical component of the 
velocity of the center of gravity of the ejection. During the 
upward moving phase, the center of gravity is decelerated 
nearly constantly, until its initial speed of Vz^Tbkm/s at 
t=40min has decreased to zero at t=46min, reaching the 
apex of the trajectory at a height of 21.5 Mm. The ejection 
then continues (roughly) along the magnetic field lines and 
falls down with a roughly constant acceleration, reaching 
about 55km/s at t=51min, just before it hits the ambient 
transition region. 

The deceleration and acceleration of the ejection are de- 
rived from a linear fit to the velocities in the time-periods 
t=40min to 45 min and t=46.5min to 50.5 min, respec- 
tively, during the rising and descending phases (Fig. [71 left 
panel, gray lines). The corresponding numbers are listed in 
Table [1 

After the initial explosion-like acceleration, the ejection 
basically follows a ballistic flight. During the rising phase, 
the front of the ejection is stronger decelerated than the 
center of gravity (^255 m/s^ for t=39.5-44 min, not shown 
in Fig. [7]) because the ejection collects the material from the 
ambient corona at its front. Consequently, the shape of the 
ejection turns from tailed to ellipsoidal. 

During the descending phase, the acceleration of the 
ejection is well below the solar surface gravitational acceler- 
ation (Table [2|). This is consistent with the mass balance as 
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Table 2. Acceleration of plasma ejection (center of grav- 
ity). 



phase 


[min] 


direction 


acceleration [m/s^] 


40.0 - 


45.0 


upward 


231 


46.5 - 


50.5 


downward 


157 



Note: Derived from linear fits to the vertical velocity shown in 

Fig.m 

discussed in Sect. 15. II After the chromospheric evaporation, 
while the ejection collects the material along the trajectory, 
part of its momentum goes into the acceleration of the col- 
lected material, which is initially at rest. Consequently, the 
ejection is not accelerated as fast as a free falling body on 
its way down along the magnetic field lines. 

5.3. Interaction with the magnetic field 

The ejection roughly follows the magnetic field lines. 
However, because of its high momentum, it is capable to 
distort the coronal magnetic field (cf. Fig.[8]). This is evi- 
dent when investigating the kinetic plasma- /3 term, which 
can be defined as the ratio of the kinetic energy density 
and the magnetic energy density i-e., 

2 

pv 

/^kin = M0-g2" • (1) 

Here, /io is the magnetic vacuum permeability, p and v are 
the mass density and velocity of the plasma ejection (at its 
center of gravity), and B is the magnetic field strength. 

In Fig. [71 (right panel) /3kin at the center of gravity of 
the ejection is shown as a function of time. Evidently, /3kin 
is approximately unity and can therefore displace the mag- 
netic field, i.e., induce electric currents. These currents are 
visible as a tail behind the ejection in the current maps 
(right panels of Fig. [1] and lower left panel in the movie at- 
tached to Fig. lA.ip . This means that the plasma ejection 
bends the magnetic field. 

Since the non-ideal MHD equations are solved in the 
numerical experiment, i.e., diffusivity is taken into ac- 
count, the ejection will start moving through the magnetic 
field lines. The characteristic time for resistive diffusion is 
T = ^^/t], where £ is the length scale over which the ejection 
moves across the magnetic field, and r] is the magnetic diffu- 
sivity. Consequently, i = (r^r)^/^. Thus, for a typical flight 
time r = 5 min and a magnetic diffusivity r] = 10^^ m^/s 
used in the numerical experiment, our estimate implies that 
the trajectory of the plasma ejection can be shifted perpen- 
dicular to the direction of the magnetic field by ^ ~ 1 ... 2 
Mm. This number is consistent with the offset of the ejec- 
tion trajectory perpendicular to the magnetic field lines (cf. 
Fig.©. 

6. Visibility in AIA/SDO 

As outlined in the introduction, cool structures are ubiqui- 
tously present in the hot corona. So far, we have focused 
on emission in O VI and Mg x to represent the hot and 
cool plasma. The Atmospheric Imaging Assembly (AIA; 
iLemen et al.L l201l[ ) on the Solar Dynamics Observatory 
(SDO) is an appropriate instrument to observe ejections 
as those investigated in the present paper. The channels 



30 p 

S5 - 




Fig. 8. Magnetic field configuration at timestep t=15min 
(well before the plasma ejection lifts off) together with tra- 
jectory of the plasma ejection (center-of-gravity) overplot- 
ted as a red line. The map at the bottom shows the vertical 
magnetic field at the bottom boundary of the simulation 
box scaled from -2000 G (blue) to +2000 G (red). 



centered at 304 A and 171 A are dominated by lines of He ll 
and Feix/x, respectively. Their temperature contribution 
functions peak at logT [K] = 4.9 and 5.9, i.e., not very far 
from the Ovi, and close to the Mgx formation tempera- 
ture (cf. Table [1]). Consequently, the images synthesized for 
the AIA channels displayed in Fig. [9] do not show significant 
differences from those in O VI and Mgx presented in Figs.[T] 
and ED 

To synthesize maps as AIA would observe the plasma 
ejection, we used the temperature response functions for the 
AIA channels as provided in Solar Soft These are based 
on synthesized emission line spectra folded with the wave- 
length response of the filters. Applying the temperature re- 
sponse function at each grid point, multiplying by the den- 
sity squared and then integrating along the line of sight pro- 
vides the count rate to be expected in the respective band, 
i.e., the digital number DN/s/pixel. In Fig.[9]these maps are 
displayed for horizontal integration (along y) for the 304 A 
band that shows the cool plasma at transition region tem- 
peratures and for the 171 A band that displays the hotter 
plasma at coronal temperatures. The images are degraded 
to match the spatial sampling of AIA (0.6 arcsec/pixel cor- 
responding to 435 km/pixel). The dynamic range of the 
AIA synthetic images in Fig. [9] is 1000 (304 A) and 300 
(171 A), which is consistent with the estimates for the AIA 
dy namic range as p rovided in the AIA instrument paper 
(Le men et al.L[2QTll ). Also, the absolute values of the count 
rates of approximately 100 DN/s/pixel are consistent with 
real solar observations of AIA. 

The plasma ejection is hardly visible in the 171 A band 
(right panel of Fig. [9]), it basically appears as a dark bubble 
with a bright ring. The reason for this is that the ejection is 
cool in the middle with the temperature rising to its circum- 
ference. However, we cannot really expect to see this ejec- 
tion in the 171 A channel in real observations. In the present 
case, the line-of-sight integration is over 50 Mm only, i.e., 
the length of the computational box. On the real Sun, struc- 

^ http:/ /www. lmsal.com/solar soft/ 
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AIA304Aband - DN s"^ pix"^ AIAlVlAband - DN s"^ pix"^ 

1 10 100 1000 10 100 1000 




X [Mm] X [Mm] 

Fig. 9. Synthesized images as they would be observed at the limb by AIA/SDO in the 304 A and 171 A bands dominated 
by the Hell line and the Feix/x lines, respectively. These snapshots shows the ejection at t=47min, i.e., at the same 
time as the middle row of Fig.[T] and as Fig. lA.ll The original synthesized images are downgraded to match the spatial 
resolution of AIA/SDO. The intensity is given in expected AIA counts per pixel and second. 



tures over a far greater length would contribute along the 
line of sight. Thus, the weak signature of the plasma ejec- 
tion can be expected to be drowned in the background. 
If the density of the ejection is sufficiently high and if the 
temperature is sufficiently low, it might significantly absorb 
light at the wavele ngths of the 171 A channel through free- 
bound absorption (|Anzer fc Heinzell , l2005f ). However, even 
then, ample structures in the foreground might contribute 
to the emission along the line of sight on the real Sun and 
cloud the visibility of the ejection in the 171 A band and 
the other coronal channels of AIA. 

The situation is different for the 304 A channel. This 
band mainly displays cooler material at less than 100.000 K 
and shows a clear signature of the plasma ejection during 
its flight. To our knowledge, no reports have been published 
so far on similar structures on ballistic orbits through the 
corona. However, it must be noted that SDO is fully opera- 
tional since last summer only. We therefore propose to look 
for these structures in the AIA 304 A data to confirm the 
co unterparts of th e (even cooler) structures seen in Ly-a 
bv lSchriiverl (|2QQl[ ). which were discussed in the introduc- 
tion, and which are consistent with the plasma ejection de- 
scribed here. When analyzing the AIA 304 A data one will 
have to struggle with the confusion of structures through 
the integration along the line of sight. Therefore, we can- 
not estimate yet if the mechanism proposed in this paper 
will be applicable for a wider range of processes seen in the 
304 A AIA channel. However, we hope that the very high 
temporal resolution of AIA will allow us to see these cool 
ejecta flying in the background through the corona. 

7. Conclusions 

We have analyzed a numerical experiment that describes 
the solar atmosphere above a small active region in which 
— without external triggering — a bubble of dense cool 
plasma is accelerated upward into the corona, moves along 
the magnetic field lines, and eventually falls back down onto 
the surface at a different place. The ejection is a direct con- 



sequence of the increased heating just above the chromo- 
sphere (and in the low transition region) for a short period 
of time. This can be considered an explosion, where mag- 
netic energy is transformed effectively to thermal energy. 
As a consequence of the explosion, a pressure perturba- 
tion travels in a ballistic fashion away from the heating site 
along the magnetic field lines. 

While the heating rate in our model is based on mag- 
netic processes (i.e., braiding of field lines and Ohmic dissi- 
pation of the induced currents), the transient event itself is 
basically a hydrodynamic phenomenon. The situation here 
is different from reconnect ion events, where plasmoids are 
driven away from the reconnection site by direct magnetic 
forcing. For the case described in this study, the increased 
heating leads to an explosion- type event, and the pressure 
gradient pushes the material up into the corona. 

The increased heating rate along the magnetic field lines 
that eventually leads to the ejection also causes (a more 
gradual) evaporation of chromospheric material before and 
during the event. Thus, the bubble that flies through the 
corona and collects mass on its way contains plasma evap- 
orated from both footpoints of its trajectory. 

The ejection has sufficient kinetic energy to deform the 
magnetic field on its fiight through the corona. This leads to 
current sheets in front of the ejection that propagate faster 
than the ejection and lead to heating (and evaporation) of 
plasma ahead of the ejection. 

Because the ejection consists of cool plasma, it can- 
not be expected to be observed in coronal images of the 
Sun that display plasma at temperatures around and above 
10^ K, and which are available through the 171 A channel 
of SDO/ AIA. The ejection might be visible in data that 
image lower temperature plasma, such as the 304 A AIA 
channel. However, more detailed studies will be needed to 
ascertain if these events are visible considering the confu- 
sion of structures along the line of sight. 

The plasma ejection observed in our numerical experi- 
ment is consistent with the type (3) brightenings as seen in 
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the TRACE Ly-a channel bv lSchriiverl (|200l[ ) and discussed 
in the introduction. Even though the plasma seen in Ly-a 
is most probably cooler than the ejection described in this 
work, our proposal of an ejection driven by an explosion- 
type event might be a good explanation for structures flying 
on ballistic orbits through the corona. Further comparisons 
between numerical experiments and observations will have 
to show to what extent this is an ubiquitous mechanism 
for transient events seen through cool plasma in the hot 
corona. 
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Appendix A: Temporal evolution of the plasma 
ejection 

The temporal evolution of the plasma ejection is best dis- 
played by a movie showing the emission as well as the 
heating rate integrated through the computational box. To 
mimic the situation at the limb, a horizontal line of sight is 
chosen, e.g., the ^-direction. The trajectory of the ejection is 
tilted by approximately 20° with respect to the x-direction 
(cf. Fig. [8]), thus, these projections nicely show the ejection 
flying through the corona. 

The intensities I{x^z]t) = e{x^y^ z]t) dy^ are dis- 
played in the top row of Fig. lA.ll where e is the emissivity 
at a given grid point x^y^z and time t. To avoid Moire pat- 
terns when displaying the images, the values for density and 
temperature were interpolated on a grid with higher spa- 
tial resolution before the emissivities were calculated using 
CHIA NTI (see als o Sect. 3]). T his procedure is similar to 
that of lPeter et all (|2QQ4 l2QQ6l ). The intensities in Fig.ET] 



are displayed on a logarithmic scale with dynamic ranges 
of 3000 (O Vl) and 50 (Mgx) to achieve a better contrast. 

The bottom panels in Fig. lA.ll and the movie show the 
current j, a measure of the Ohmic heating rate (cxj^), in 
two different normalizations. The lower left panel shows p 
normalized by the average trend with height, i.e., the line- 
of-sight-integrated currents {Jy • • • ^v) averaged along the 
x-direction and in time, (. . .)x^t^ 

J y, z; t) dy 

This normalization is necessary to see structures in 
the heating rate (ocj^), which when averaged horizontally, 
drops roughly exponentially with a scale height of about 
5Mm (|Bingert Peteil l201lD . 

To emphasize the temporal variability, the heating rate 
was normalized at each grid point in the 3D computational 
domain by the temporal average, (. . .)^, at the respective 
grid point. This normalized quantity is then integrated 
along the line of sight, i.e., 

jI {x,z;t) = [ . d^ . (A.2) 

Jy {f{x,y,z;t) )t 

The term basically shows fluctuations along the line of 
sight, which is why it is indexed F. In the lower right panel 
of Fig. lA.ll this quantity is plotted using a rainbow color 
table, where violet represents low values of and red rep- 
resents high values. 

A potential magnetic fleld is used as an upper bound- 
ary condition above the computational domain in the 3D 
MHD model. Of course, the magnetic fleld inside the box 
is computed self-consistently based on the MHD equations. 
As a result, increased currents are observed near the top 
boundary. These currents have no impact on the magnetic 
structure and plasma heating in the main part of the box. 
They are small compared to those at lower heights, which 
are responsible for heating the coronal plasma (note that 
the average currents decrease exponentially with height). 
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Fig. A.l. Side view of the 3D computational box showing the synthesized emission in Ovi and Mgx (top row), the 
normalized currents jj^ and the current fluctuations jp (bottom row; see Eqs. lA. II and IA.2I and Appendix [Al for details). 
Each panel covers 50 Mm in the horizontal and 30 Mm in the vertical direction. This snapshot shows the ejection at 
t=47min, i.e., at the same time as the middle row of Fig. [T] and as Fig.[9l The arrow points at the location of the high 
energy input that causes the ejection. A movie showing the temporal evolutio n over 1 h is av ailable in the online edition. 
The movie is also available at |http: / / www.mps.mpg.de/data/outgoing/ peter / papers / plasma-ejection/blob. mpg (6 MB). 



